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Ultrasonic attenuation was measured for polycrystalllne samples of nickel 
and copper with various grain-size distributions produced by heat treatment. 
Attenuation as a function of frequency was determined for a sample having a 
known mean grain diameter. Once this function was determined, It could be 
scaled to determine the mean grain size of other samples of the same material 
with different mean grain diameters. These results were obtained by using 
broadband pulse-echo ultrasound In the 25- to 100-MHz frequency range. The 
results suggest an ultrasonic, nondestructive approach for verifying heat 
treatment of metals. 


INTRODUCTION 

An extensive effort has been made to relate the observed ultrasonic atten- 
uation to the mean grain diameter In polycrystalllne solids (refs. 1 to 7). 
Although the experimental and theoretical results agree for some polycrystal- 
llne materials, many materials are so complex that they are not adequately 
described by existing theory (refs. 8 to 13). Because structural materials 
generally exhibit complex microstructures, verification of their microstruc- 
tures Is of high Interest. As a practical matter, determining the mean grain 
size of both the simple and complex systems by nondestructive ultrasonic atten- 
uation Is preferable to the usual time-consuming, expensive, and Inherently 
destructive metallographlc methods. 

Most polycrystalllne solids are a collection of densely packed, Inhomoge- 
neous, and morphologically complicated scatters. Previous Investigators have 
usually assumed that a particular mode of scatter attenuation (Rayleigh, 
stochastic, etc.) Is dominant when comparing experimental and theoretical 
results. However, In many complex systems both Rayleigh and stochastic scat- 
tering may occur for a given range of frequencies (or wavelengths). 

In general Rayleigh, stochastic, and diffusive scattering may occur simul- 
taneously. The relative contribution of each scattering mechanism to the total 
observed attenuation Is unknown for a complex microstructure. In this study no 
attempt was made to attribute attenuation to a particular scattering mode. We 
merely observed that a significant factor In comparing changes In microstruc- 
ture Is a ratio containing wavelength and mean grain diameter. It will be 
shown that any arbitrary wavelength can be used to establish a mean-grain-si ze- 
related scale factor between two polycrystalllne microstructures. 
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THEORY 


Expressions for ultrasonic attenuation coefficient a as a function of 
frequency f have been developed for three regimes: the Rayleigh scattering 

regime, where wavelength x greatly exceeds mean grain size D (X » D) 

= C R D 3 f 4 , Np/cm (1) 

the stochastic scattering regime, where X ~ 0, 

a S = c s ^ f2, N P /cm ( 2 ) 


and the diffusion scattering regime, where X « D, 

C . 

a, ■ — , Np/cm (3) 

d D 


Here, Cr, C s , and are constants that depend on material parameters such 
as elastic constants and longitudinal and shear velocities. Typical dimen- 
sions for a are nepers per centimeter. Equations (1) to (3) are generally 
valid only for homogeneous, regular, and equlaxed systems (ref. 10). 

Note that ea£h of the expressions (eqs. (1) to (3)) contains a coefficient 
that depends on D, the mean grain slze^ Assume that the constants Cr, C s , 
and C,j are unaffected by changes In 0. Also assume that the shape of the 
grain size distribution function and grain morphology are unaltered by changes 
In D. Given these assumptions, the microstructure would simply change by a 
scale factor corresponding to the Increase (or decrease) In T). Substituting, 
V/x, where V denotes velocity, for f and multiplying equations (1) to (3) 
by wavelength x yields a new set of equations In which D/x Is the key vari- 
able. Thus, given that velocity Is not a function of frequency, we have for 
Rayleigh scattering 


a 


X 


s 


I 




3 

, Np/wavelength 


for stochastic scattering 



Np/wavelength 


for diffusion scattering, 


a X = C d (~^J’ N P/ wavelen 9 th 


(4) 


(5) 


( 6 ) 


where Cr, Cj, and C<j are constants. The attenuation coefficient Is now 
expressed In nepers per wavelength for each regime. 
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Evaluating equation (5) at some arbitrary wavelength \ 0 for a system 
So having a mean grain diameter Dq (fig. 1(a)), we obtain the attenuation 
coefficient 


(v 


^0’N)' C « 


= a g, Np/wavelength 


( 7 ) 


If the mean grain size and the wavelength are doubled (fig. 1(b)), we obtain 
from equation (5), 


r' flE 

“k ■ C s U, 


2Do»2^o'**i 


otg, Np/wavelength 


( 8 ) 


or the same numerical value for the attenuation coefficient as obtained In the 
previous case (similar results can be obtained by using eq. (4) or (6)). This 
Is shown graphically In figure 2(a). By varying the arbitrary wavelength \ 0 , 
we can construct the attenuatlon-versus-wavelength functions So and S 
(fig. 2(b)). Because the mean grain size of system S Is also arbitrary, we 
can construct the attenuatlon-versus-wavelength function for the N th system 
S n having the mean grain size D n . 

A scale factor can be defined as 


a = = — = N** 1 scale factor (9) 

n \) B 0 

where O n /Do Is the ratio of the mean grain size of the N th system to the 
mean grain size of the initial system. The scale factor Is also represented 
as the ratio of wavelengths for which constant attenuation <*o Is observed 
between two systems (l.e., the ratio of wavelength \ n of the N th system 
to the wavelength \ 0 of the Initial system at constant attenuation). 

The curve Sq In figure 2(b) Is the generating function from which the 
other curves can be obtained. To construct the generating function, a fourth- 
degree polynomial In frequency was fitted to the attenuation data by using the 
method developed by Marquardt (refs. 14 and 15). Here the attenuation Is 
expressed as 


* 0 ( f ) 



Np/cm 


( 10 ) 


where A m are constants. The terms In the expansion shown In equation (10) 
are not directly related to equations (1) to (3). Any other analytical expres- 
sion that fits the attenuation data would be equally valid. The preceding 
curve-fitting function Is made a function of wavelength by combining 
equation 10 with 
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V = fx 


( 11 ) 


where V is the velocity of sound, to obtain 


v*-> 


E 

m=l 




m 


Np/cm 


( 12 ) 


where are constants. By multiplying equation (12) by the variable wave- 
length, we obtain the generating function 



Np/wavelength 


(13) 


The theoretical attenuation for other systems (of larger grain size) Is 
obtained from their respective generating functions by scaling the wavelength 
as Indicated In figure 2. The scale factor cr Is varied until the generating 
function for the Initial system with the scaled wavelengths reconstructs the 
attenuation for the N th system. If the mean grain size of the Initial system 
Is known, the mean grain size D n for the N th system Is given by 
equation (9). 


EXPERIMENT 
Material Samples 

Two materials, nickel 200 and copper, were used In this Investigation. 

The nickel 200 was 99.5 percent pure and the copper was 99.99 percent pure. 

Both were obtained as commercially available rods. The nickel 200 was then 
cold rolled 50 percent and annealed In air at 690 + 9 °C and cold rolled again 
50 percent. The copper sample was cold rolled 60 percent and annealed at 
329 °C In argon for 5 min. Three nickel samples, cut from the rolled stock, 
had typical dimensions of 1 by 1 by 0.2540 cm. Five copper samples, cut from 
the rolled stock, had typical dimensions of 1.2 by 0.8 by 0.18 cm. Each sample 
was polished to a mirror finish with a final polish of 1 .0-nm-dlameter diamond 
slurry. 


Microstructure Modification 

The mean grain size was Increased In both the nickel and copper samples 
by annealing (refs. 16 and 17). Three nickel samples were annealed In air at 
800 °C for 1, 2, and 5 min, respectively. Four of the five copper samples 
were annealed In argon for 5 min at 466, 510, 700, and 800 ®C, respectively. 
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Attenuation Measurement 


Before modifying the Initial grain size, we made attenuation measurements 
on each sample. The methodology for making precision attenuation measurements 
Is described In reference 18. These Initial attenuation measurements were made 
to verify that all of the nickel and copper samples were nearly acoustically 
Identical to the other samples of the same material. At 50 MHz there was 4 and 
6 percent variation In the attenuation measurement among nickel and copper 
samples, respectively. These variations Included the effect of the Imprecision 
In attenuation measurements that Include the frequency-dependent reflection 
coefficient In the calculation (ref. 18). 

After the mean grain size had been modified, attenuation measurements were 
again made on each sample. Since the attenuation coefficient (at fixed fre- 
quency) generally Increases with mean grain size, samples exhibiting extremely 
high attenuation were ground thinner so that the uncertainty In the attenuation 
measurement was minimized (ref. 18). The results were reproducible with typi- 
cal variations In measurement similar to those for the Initial attenuation 
measurement. 


RESULTS 

Experimental Generating Functions 

In general, for a fixed frequency, the attenuation coefficient Increased 
with annealing time (fig. 3). The solid curve Is the best-fit curve (l.e., 
generating function) obtained by using equation (10), for the data for the 
nickel sample annealed for 1 min. The data shown In figure 3 reappear In 

figure 4 as a function of wavelength. The solid curves for the samples 

annealed for 2 and 5 min were determined from the generating function by scal- 
ing the wavelength, as discussed earlier. The observed and theoretical results 
are shown In figure 5 as a function of frequency. 

For fixed frequency the attenuation coefficient Increased with annealing 
temperature (fig. 6) for the copper samples. The solid curve for the copper 
sample annealed at 329 °C Is the generating function for these samples. The 
data shown In figure 6 reappear In figure 7 as a function of wavelength. The 
solid curves were determined from the generating function by scaling the wave- 
length. The observed and theoretical results are shown In figure 8 as a func- 
tion of frequency. 

The scale factors for the data shown In figures 4 to 8 were determined by 

scaling the generating function to fit the frequency regime for which the most 

accurate data were obtained. This regime, although well defined (ref. 18), was 
different for each sample and depended on the sample thickness, the attenua- 
tion, and the frequency-dependent reflection coefficient. 

Deviations existed between the experimental and theoretical results 
(figs. 4 to 8) at both the low and high frequencies (l.e., long and short wave- 
lengths). These deviations were due to the greater uncertainty In the measure- 
ment of attenuation at these extreme frequencies (ref. 18). 
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Mean Grain Size Ratios 


Several methods are available for determining the mean grain size 
(refs. 12 and 19). For regular systems the standard line-intercept method 
appears to be adequate. However, this method Is not applicable to systems with 
structurally complicated morphologies. More advanced theories that are based 
on spherical grains also are Inherently not applicable to real, complicated 
systems. An alternative technique, the comparison method (ref. 19), Is free of 
the theoretical bias Imposed by the above-mentioned theories. 

The simple comparison method allows for the grain structure to be visually 
compared with an ASTM standard chart. This technique Is not directly applic- 
able to complicated systems ( 1 . e . , the standard ASTM charts are for regular 
structures). We can modify the comparison method so that we may compare 
systems with dimensionally arbitrary scales but of similar complexity. This Is 
done by taking photomicrographs of similar systems of different scales and sub- 
sequently magnifying the photomicrograph for the smaller system by the appro- 
priate scale factor between the two systems. The magnified photomicrograph and 
the photomicrograph for the system of larger scale should appear similar. This 
Is shown graphically In figures 9 and 10 for the nickel and copper samples, 
respectively. (Ultrasonic attenuation was measured along the vertical axis.) 

From figure 4, the scale factor between the photomicrographs for the 
nickel systems annealed for 2 and 5 min (figs. 9(a) and ( b) ) Is 1.85. 

Figure 9(a) Is shown magnified 1.85 times In figure 9(c), which Is dimension- 
ally and morphologically similar to figure 9(b). Similar results are shown In 
figure 10 for the copper samples annealed at 510 and 800 °C. Comparisons with 
the other systems show strikingly similar results. 


DISCUSSION 

The success of scaling a generating function over a wide frequency range 
to other systems of similar complexity but of different scale was quite remark- 
able In that only one free parameter was needed. This parameter Is directly 
related to the relative scale between systems. Ideally, one would like to use 
(scale) a theoretical expression for ultrasonic attenuation that Includes 
multiple grain-boundary scattering for morphologically complicated systems. 
However, such an expression Is currently unavailable since the theory of grain- 
boundary scattering for complicated systems remains underdeveloped (ref. 10). 

The lack of a comprehensive theory can be bypassed by measuring the atten- 
uation for an Initial system as a function of frequency. The effects of grain 
distribution, morphology, orientation, and grain Interaction are Implicit. 

This experimentally measured attenuation also Includes attenuation due to dis- 
location damping (refs. 6 and 20), thermoelastic effects (ref. 20), etc. These 
nonbarrier attenuation mechanisms will not scale as the grain-boundary scatter- 
ing terms (eqs. (4) to (6)) do. We have assumed In the preceding analysis that 
these unscalable attenuation mechanisms are negligible relative to the observed 
attenuation. This, however, will not be true for materials that are In a 
mechanically worked state, where the dislocation densities are expected to be 
large (ref. 16). 
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A major requirement for the scaling method to be successful Is that the 
systems of Interest must have similar morphology. This Is a somewhat subjec- 
tive Judgment not quantitatively described. The lack of a comprehensive theory 
for determining the mean grain size for a real, complicated system has left us 
with no alternative, other than the comparison method, for making quantitative 
comparisons between theory and experiment. Although magnifying the photomicro- 
graphs (figs. 9 and 10) served to make these comparisons, It Is not the best 
method for comparing grain structures. A more precise method would be to use 
a continuously variable magnification metallograph. 

The theories for both ultrasonic boundary scattering and grain size dis- 
tribution are Insufficiently developed to provide quantitative results for com- 
paring experiment and theory. This problem Is temporarily resolved by using 
the observed attenuation and the observed microstructure to show this agree- 
ment. The full comprehensive theory for the ultrasonic attenuation that 
Includes grain-boundary scattering for morphologically complicated systems and 
allows for Interacting grains Is expected to scale as the experimental results 
Indicate here. 

The simple concept of mean grain size Is likely to be Inappropriate for 
both acoustic and microstructure theory. A more complete description of the 
microstructure morphology would Include the distribution of grain size 
(ref. 12) and grain boundary topology. This more complete description should 
be Included In determining of the microstructure morphology and ultrasonic 
attenuation theory. Nevertheless, even In the form described In this paper 
the approach for determining microstructure scale Is likely to be useful In 
verifying the heat treatment of most polycrystalllne materials. 


CONCLUSIONS 

It has been shown that the ratio between wavelength and mean grain 
diameter Is a significant factor In comparing changes In microstructure. This 
ratio can be used to determine a scale factor between two polycrystalllne 
materials with similar morphologies but different mean grain diameters. These 
results were achieved by determining the attenuation as a function of frequency 
without reference to any particular conventionally accepted scattering mode. 

By this approach It was unnecessary to associate a particular wavelength with 
mean grain diameter. This allowed a generalized formulation for scaling mean- 
grain-diameter variations among samples of a polycrystalllne material that had 
undergone a range of heat treatments. 
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ATTENUATION COEFFICIENT, 
alf), Np/cm 



WAVELENGTH, \ cm 

<b) Systems S 0 , S, and S N . Curve labeled S 0 
is generating function. 


Figure 2. - Attenuation coefficient as function 
of wavelength. 
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Figure 3. - Attenuation coefficient as function of 
frequency for nickel samples. Annealing temper- 
ature, 800 °C. Solid curve is generating function. 
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Figure 4. - Attenuation coeffi- 
cient as function of wave- 
length for nickel samples. 
Annealing temperature, 

800 °C. 
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Figure 5. - Comparison between observed and theo- 
retical attenuation coefficients as function of fre- 
quency for nickel samples. Annealing temperature, 
800 °C. 



Figure 6. - Attenuation coefficients as function 
of frequency for copper samples. Annealing 
time, 5 min. Solid curve is generating func- 
tion. 
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Figure 7. - Attenuation coefficient as function of wavelength for copper 
samples. Annealing time, 5 min. 
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Figure 8, - Comparison between observed and theoretical atten- 
uation coefficients as function of frequency for copper sam- 
ples, Annealing time, 5 min. 
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SCALE FACTOR, 1.85 


100 pm 

(a) Annealed for 2 min at 800 °C. 


(b) Annealed for 5 min at 800 °C, 


(c) 1.85 Magnification of figure 9(a). 


Figure 9. - Photomicrographs of nickel sample. 


SCALE FACTOR, 1.75 


50 pm 

(a) Annealed for 5 min at 510 °C. 


(b) Annealed for 5 min at 800 °C 


(c) 1. 75 Magnification of figure 10(a). 


Figure 10. - Photomicrographs of copper samples. 
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